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UNSTEADY LIFT FORCES GENERATED BY VORTEX SHEDDING ABOUT 


A LARGE, STATIONARY, AND OSCILLATING CYLINDER AT 
HIGH REYNOLDS NUMBERS 
By George W. Jones, Jr. 

NASA Langley Research Center 

ABSTRACT 

A wind-tunnel study has been made of the unsteady lift (crosswind) forces 
on a large, two-dimensional cylinder with axis perpendicular to the flow. The 
forces were measured at Reynolds numbers from 0.J6 million to 18.7 million at 
Mach numbers up to 0.6 with the cylinder held both stationary and oscillated 
laterally. The mean-drag and the Strouhal number of the unsteady lift forces 
on the stationary cylinder followed the trends established by previous inves- 
tigators. At higher, hitherto -unexplored Reynolds numbers, the mean-drag was 
approximately 0.54, and the Strouhal number was approximately 0*3* The 
unsteady root -mean- square lift force on the stationary cylinder was defined 
over the flow parameter range investigated. For the oscillated cylinder, an 
unsteady lift force due to motion was found to exist at cylinder frequencies 
near the vortex- shedding frequency of the stationary cylinder. This unsteady 
lift due to motion was found to be a destabilizing aerodynamic damping force 
at frequencies below the stationary cylinder vortex-shedding frequency which 
shifted abruptly to a stabilizing damping force at frequencies above the 
vortex- shedding frequency. 



INTRODUCTION 


Recently, the response to ground winds of large flexible structures such 
as metal smokestacks and cylindrical space vehicles standing on the launch pad 
has led to considerable interest in the spontaneously generated unsteady aero- 
dynamic forces on circular cylinders at high Reynolds numbers. These unsteady 
forces are produced by vortices shed alternately from the sides of the struc- 
ture and, for lightly damped structures, are capable of exciting large ampli- 
tude vibratory response. The unsteady forces associated with vortex shedding 
are not well understood or amenable to analysis. Many investigators have made 
experimental studies such as given in references 1 to 11 which have given some 
understanding of the mechanism of vortex shedding, particularly at lower 
Reynolds numbers. However, at higher Reynolds numbers, there is little inform- 
ation on the nature of these aerodynamic forces. In particular, the question 
of how these forces are affected by the motion of the structure as it deflects 
in response to the aerodynamic loads has needed investigation. 

A wind-tunnel investigation has, therefore, been conducted on a large 
circular cylinder in a two-dimensional flow at Reynolds numbers from O.36 mil- 
lion to 18.7 million. The cylinder was instrumented to read directly the mean- 
drag and unsteady lift forces. In addition to being fixed, the cylinder could 
be oscillated over a range of frequencies and amplitudes. This oscillation 
capability was used to investigate the effects of cylinder motion on the aero- 
dynamic forces generated. Figure 1 presents a comparison of the Reynolds 
number and nondimens ional frequency- of- oscillation ranges of the present inves- 
tigation with those of previous investigations. This figure shows the present 
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investigation has greatly extended the range of Reynolds number and nondimen- 
sional frequency of oscillation over -which the aerodynamic forces on cylinders 
have been studied. 

The investigation was a joint effort by NASA., Marshall Space Flight Center, 
and NASA, Langley Research Center, together with The Martin Company, Baltimore. 
The Martin Company, under contract, supplied the model and participated in the 
tests and some data reduction. Results of this investigation, based on pre- 
liminary data reduction, were given in reference 12. The present paper will 
first briefly review some of the major results obtained in the study and show 
some comparison with other investigators which was lacking in reference 12. 

Then the paper will concentrate oh the effects of cylinder motion on the 
unsteady lift forces. It should be pointed out that the data on cylinder- 
motion effects presented in reference 12 indicated correct trends, but were 
.quantitatively in error. 

SYMBOLS 

A projected area of instrumented section of cylinder 

Cj) mean-drag coefficient for instrumented section of cylinder 

G L,Imag imaginary component (in phase with cylinder velocity) of the 

aerodynamic-lift-coefficient vector at the frequency of 
oscillation of the cylinder 

0^ r 6 real component (in phase with cylinder displacement) of the 

aerodynamic-lifb -coefficient vector at the frequency of 
oscillation of the cylinder 

Cl rillR root -mean- square value of unsteady-lift coefficient on 

f instrumented section of cylinder 

C-p base average of section-pressure coefficients at rear of cylinder 

’ between cp = 150° and cp = 210° 

D diameter of test cylinder, 3 feet (O.^HM- meter) 
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frequency, Hz 

forcing frequency of cylinder oscillation, Hz 
vortex-shedding frequency, Hz 


instantaneous amplitude of cylinder displacement under forced 
oscillation, h = £lq cos oot 


Serut on-type aerodynamic derivative of real component of 
unsteady-lift-force-coefficient vector at frequency of 


cylinder oscillation. 


~^L,Re 

2 S 2(§) 


amplitude ( 0-peak) of cylinder displacement 


Scruton-type aerodynamic derivative of imaginary component 
of unsteady-lift-force-coefficient vector at frequency of 


cylinder oscillation, 


" c L,Imag 



Mach number 

free-stream dynamic pressure, — pV 2 

Reynolds number based on cylinder cross-section diameter 


aerodynamic Strouhal number, 


££ 

V 


time 

velocity of test medium 
viscous damping ratio 
density of test medium 

angular location of static pressure orifice ■with respect 
to direction of free-stream flow 

power spectral density of unsteady-lift- force coefficient as 

00 

$C L (S)dS 


function of Strouhal number, (Lift) 2 = (qA) 2 


0 


circular frequency 



APPARATUS ARD TESTS 


The Langley transonic dynamics wind, tunnel used for the investigation 
has a test section approximately 1 6 feet (4.877 meters ) square with cropped 
corners and the unique ability to use either air or Freon 12 as a test medium. 

A photograph of the cylinder model mounted in the tunnel is shown in figure 2. 
Although the model cross-section area appears large in relation to the tunnel 
cross-section area, the estimated solid-blockage interference was negligible, 
and no wall- interference corrections have been made to the data. 

The cylinder model has a 3-foot (0.914— meter.) diameter and spanned the 
tunnel with its longitudinal axis perpendicular to the air flow. The model 
structure consisted of an inner support cylinder, an outer cylinder divided 
into sections, a variable-stiffness suspension system, and a hydraulic shaker 
system. The inner cylinder, a seamless aluminum tube, carried the load and 
extended through the tunnel floor and ceiling to attach to the suspension 
system. The outer-cylinder sections were of lightweight aluminum honeycomb 
sandwiched between thin aluminum sheets. The top and bottom outer-cylinder 
sections were rigidly fastened to the inner cylinder with spacers. A 7-foot 
(2.1336-meter) center section of the outer cylinder was attached to the inner 
cylinder by force transducers, oriented so as to read directly the integrated 
forces on the center section in the direction of and at right angles to the 
flow. The 0.02- inch (O.o4o8-cm) gaps between the fixed top and bottom sections 
of the outer cylinder and the instrumented center section were originally 
designed with a labyrinth seal which proved unsatisfactory, and was supple- 
mented by filling the gaps with a silicone rubber compound. This flexible 
rubber seal allowed relative motion between adjacent cylinder sections and 
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transmitted in shear only about 1 percent of the load on the center section. 
The suspension system consisted of one tapered torque bar* attached to the top 
of the cylinder and an identical bar attached to the bottom of the cylinder. 
Remotely controlled clamping mechanisms could be positioned along the torque 
bars so as to change the resonant frequency of the system oyer a range from 
3 to 20 Hz. Two hydraulic shakers, one positioned at the top and the other 
at the bottom of the cylinder, were driven from a single control, and could 
either hold the cylinder fixed or impart an oscillatory driving force at right 
angles to the tunnel flow on the resonant system over the range of resonant 
frequencies at amplitudes up to 3 inches displacement. 

A Martin-developed instrumentation system called the Inertia Compensated 
Balance was used to measure directly the aerodynamic forces on the 7-foot 
(2.1336-meter) outer-cylinder section in the plane of the flow and at right 
angles to the flow. This ICB system electrically combined the signals from M 
the strain-gage force transducers with the signals from two strategically 
located accelerometers in such a manner as to cancel out the forces due to 
inertia and thus give a direct measure of the aerodynamic forces. A single 
row of static-pressure orifices measured the static -pres sure distribution 
around the cylinder at one longitudinal station, 1.1*1- diameters from the lower 
end of the model. 

During the test program, the mean- drag and unsteady lift forces were 

measured at Reynolds numbers from O.36 to 18.7 million at Mach numbers up to 

0.6 with the cylinder held fixed and with the cylinder oscillated at model 

f h D . . 

reduced frequency, — , from 0.06 to 0.50. Static pressures were measured 
ondy with the cylinder held fixed. 
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SELECTED RESULTS AND DISCUSSION 


As stated in the introduction, the results presented here will briefly 
review some of the mean-drag and unsteady lift data covered in reference 12, 
giving comparisons with other data not presented in reference 12. Following 
this review, the results for the aerodynamic forces on the moving cylinder will 
be presented since these results, as presented in reference 12, are incomplete 
and quantitatively in error. 


Mean-Drag Data 

The mean-drag data measured on the cylinder at M $ 0.2 are presented 
in figure 3 and compared with data from previous investigators. These data 
are particularly meaningful for ground winds on launch vehicles or similar 
large cylindrical structures where Mach numbers are low but Reynolds numbers 
are high. The trends from four different sources agree well; however, at given 
Reynolds numbers the differences in measured drag coefficients are appreciable. 
The greater drag on Roshko f s sandblasted pipe cylinder, as compared with the 
drag on the smooth cylinder of this investigation, was probably caused by the 
difference in roughness of the two cylinders. The data show that, for a smooth 
cylinder, the drag coefficient increases from 0.15 at Rjj « 0.5 X 10^ to about 
0.54 at % ~ 3 x 10^, and then decreases slightly as Reynolds number increases 
to 10 x 10 6 . 

The variation of base pressure coefficient is presented in figure 4 for 
Mach number <0.2. Base pressure data from references 5 and 8 are also pre- 
sented for comparison. Note that these base pressures show the same trends 
as the drag coefficients of figure 3* 



Characteristics of Unsteady Lift Measured on 
Stationary Cylinder 

Frequency content of unsteady lift .- The unsteady lift on a stationary 
cylinder can be periodic, random, or combinations of both forms. In order to 
deal with these forms, a classification of the unsteady lift on a stationary 
cylinder in terms of the frequency content has been adopted as follows j 

The usual nondimensional aerodynamic Strouhal number is defined as 
f D 

S = where f„ is the vortex- shedding frequency of the cylinder. If a 

power spectral density plot is made from a time history data sample of the 
unsteady lift force, and the frequency of the main peak in this plot is defined 
to be f s , then the frequency interval between the half-power points of this 
spectrum peak may be defined as the Strouhal bandwidth. Figure 3 shows, in 
inset form, a typical power spectral density plot with the Strouhal bandwidth 
indicated. The larger plot shows the variation of this Strouhal bandwidth 
with Reynolds number. These data are limited to M < 0.3* From study of this 
plot and the associated autocorrelation and probability -density functions 
determined from the data time history samples, the frequency content of the 
unsteady lift on the stationary cylinder may conveniently be divided into the 
three regimes indicated in figure 5* These regimes and the Reynolds number 
ranges over which they apply are wide -band random (l.l x 10^ < R^ < 3-5 x 10^), 
narrow-band random (3*5 X 10^ < Rjj < 6.0 x 10^), and quasi-periodic 
(R n > 6.0 x 10^). The wide-band random forces occur in what has been termed 
the "supercritical" Reynolds number range in references 2, 5; and 8. The 
quasi-periodic forces occur in the "trans critical" Reynolds number range of 
reference 5, and the narrow-band random forces occur in an intermediate or 
transition range between the wide-band random and quasi-periodic regimes. 
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Previous investigators have presented the frequency characteristics of the 
unsteady lift force in terms of a single aerodynamic Strouhal number rather 
than a Strouhal bandwidth. For purposes of comparison, the data of the present 
investigation are presented in this form in figure 6 . Also presented for com- 
parison are data from other investigators. The shaded areas denote ranges 
where a definite, easily defined vortex-shedding frequency occurred. In the 
Reynolds number range from about O .56 million to 2.6 million, there is no 
distinct frequency, and all investigators had difficulty in identifying a pre- 
dominant frequency. This fact may account for the considerable scatter in this 
region. It is seen that there is good agreement between the present data and 
that of Roshko (ref. 5 ) in the Reynolds number range between 5*7 million and 
8.4 million, and that in the previously unexplored region between 8.4 million 
and 18 million the Strouhal number remains relatively constant at about 0.3* 
Magnitude of unsteady lift .- The root -mean- square lift coefficients on 
the stationary cylinder at M < 0.3 and M = 0.4 are presented in figure 7 
as functions of Reynolds number. Also presented for comparison are similar 
data from Fung (ref. 2) and Schmidt (refs. 7 and 8 ). Examination of this 
figure shows that for the data at M < 0.3* the unsteady lift appears to be 
multiple valued at Reynolds numbers below about 8 million. This phenomenon 
is believed to be related to roughness asymmetries in the laminar flow regions 
on the cylinder. These asymmetries change from time to time as small particles, 
or scratches collect on the upstream surface of the cylinder during a test, and 
therefore give rise to apparent scatter in the data. Above 8 million Reynolds 
number, where the flow is turbulent, the root -mean- square lift coefficient 
measurements were repeatable and followed the single curve shown in figure 7 * 
Thus, the multiple-valued data seem to occur at Reynolds numbers corresponding 
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to those for the wide-hand and narrow-hand ranges of Strouhal numbers and merge 
to a single-valued boundary at Reynolds number where the yortex shedding is 
quasi-periodic . For Mach number 0.4, the lift coefficient measurements are 
repeatable for all points shown in figure 1, and there is a large Mach number 
effect at Reynolds numbers from 7 million to 17 million as shown by the rounded 
peak in the data at a Reynolds number of about 9 million. 

Comparative data from Fung (ref. 2) and Schmidt (refs. 7 and 8 ) are shown 
in figure 7 in the form of broad bands because of scatter in the individual 
data points. Both sets of comparative data tend to merge into the data of the 
present investigation. Fung’s data (ref. 2) approach the upper part of the 
multiple -valued region, and Schmidt’s data (refs. 7 and 8 ) merge with the 
lower part of this region. Thus, the unsteady lift coefficients obtained on 
the stationary cylinder are consistent with the trends of previous data. 

Characteristics of Unsteady Lift Due to 
Cylinder Motion 

In order to investigate the effects of cylinder motion on the unsteady 
lift force, a positive displacement sinusoidal oscillation was imparted to the 
model by the hydraulic shakers at the top and bottom of the cylinder. The 
model was driven at frequencies from 3 bo 18.7 Hz and amplitudes to 3 inches 
(0.0762 meter). The forces in the lift direction were measured on the 7-foot 
( 2 . 1336 -meter) center cylinder section using the Inertia Compensated Balance 
system, and the data thus obtained were reduced using a variety of techniques. 
The major results will be presented first in terms of the same characteristics 
used to describe the lift on the stationary cylinder. Next, by playing the 
tape recorded data samples of unsteady lift force into a special purpose analog 
computer, the in-phase and 90 ° out-of-phase (real and imaginary) components of 
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the unsteady lift force with respect to displacement were determined at the 
frequency of motion. Results from this analysis are then presented. 

Power spectral density and root-mean-square values .- Figure 8 presents 
the power spectral density of the lift force in coefficient form for two data 
points at the same Mach number and Reynolds number. For one data sample the 
cylinder was held stationary, and for the other the cylinder was oscillated at 
the aerodynamic Strouhal frequency of the stationary cylinder (S = 0.31)* A 
large increase in the peak of the spectrum for the moving cylinder is readily 
apparent. The root -mean- square lift coefficient increased from 0.0762 to 
0.1296. Thus, this comparison shows the existence of a lift due to motion. 

In order to show the ratio of this unsteady lift due to motion to the 
unsteady lift on the stationary cylinder, selected values of the root-mean- 
square coefficient of lift on the moving cylinder were normalized by corre- 
sponding values on the stationary cylinder. The results are plotted in fig- 
ure 9 as functions of nondimens ional frequency ratio and cylinder amplitude of 

motion. The frequency ratio is the nondimensional frequency of forced oscil- 
f* D 

lation, divided by the aerodynamic Strouhal number for the stationary 

cylinder. The data for this figure have been restricted to Maqh numbers less 
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than 0.3 and Reynolds numbers greater than 5*5 x 10°, where Mach number effects 
are small and the periodicity of the flow is clearly defined. 

The lift coefficient ratios of figure 9 show that when the cylinder is 
oscillated in the lift plane at frequencies far removed from the aerodynamic 
Strouhal frequency, there is no increase, or, in some eases, a slight decrease 
in lift over that of the stationary cylinder. However, for cylinder motion at 
or near the Strouhal frequency, there is a definite lift due to motion, and 
this lift increases with amplitude of cylinder motion. The maximum 
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amplification of the lift due to motion occurs when the frequency is approxi- 
mately equal to 0-99 and increases with increasing amplitude as indicated by 
the small plot at the lower right-hand side of figure 9* 

Lift force vector at frequency of cylinder oscillation .- Examination of 
data such as shown in figures 8 and 9 and of the time histories of the lift 
force on the oscillating cylinder indicated a likelihood of correlation between 
the aerodynamic forces and the cylinder motion. Therefore, an analysis was 
made using the tape recorded time histories of the unsteady lift force from 
the data samples of figure 9 and the corresponding time histories of sinusoidal 
cylinder motion. This analysis determined the aerodynamic force components 
that are in phase (real component) and 90 ° out of phase (imaginary component) 
with the sinusoidal displacement. The results of this analysis are presented 
in figure 10 which gives the real and Imaginary vector components with respect 
to the displacement of the lift force on the moving cylinder at the frequency 
of cylinder oscillation. These components have been normalized by the root- 
mean-square lift component on the stationary cylinder and are presented as 
functions of nondimensional frequency ratio and amplitude. 

The real component may be thought of as an aerodynamic stiffness parameter, 
and it behaves much as the root-mean-square unsteady lift on the oscillating 
cylinder (fig. 9) in that the real component (fig. 10(a)) of the lift vector 
is substantial only when the cylinder oscillates at or near the aerodynamic 
Strouhal frequency, and this component increases with increasing cylinder 
amplitude of motion. It should be noted that the apparent mass effects are 
not included in this real component since the Inertia Compensated Balance 
system cancels them out even though they might legitimately be included here. 
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NASA- Langley, 1967 


Figure 11. - Effect of damping on Saturn V ground wind loads. 



Similarly, the imaginary component of the lift coefficient vector 
(fig. 10(b)), being in phase with the velocity, may be regarded as ah aero- 
dynamic damping parameter. This component of the unsteady lift due to motion 
was found to be a destabilizing aerodynamic damping force at frequencies below 
the stationary cylinder vortex- shedding frequency. This force shifts abruptly 
to a stabilizing damping force at frequencies above the vortex- shedding fre- 
quency. Also, this damping component changes with amplitude of motion as 
shown in the small inset figure at the lower right. 

These real and imaginary lift coefficients can be converted into aero- 
dynamic derivatives such as used by Scruton in reference 4. The relationships 
between Scruton's derivatives and the real and imaginary lift coefficients of 
figure 10 are 


- 1 % 


c L,Re 


and 




< -'L,Imag 


For flexible structures such as a launch vehicle on its pad, a nonlinear aero- 
dynamic damping term such as -1% can give rise to a self-excited limit ampli- 
tude sinusoidal oscillation. Such oscillations have been observed on aeroelastic 
ground wind load models of launch vehicles (ref. 15 ) and on other structures 
(ref. 14). 

Figure 11 shows a typical example of such an oscillation on a Saturn V 
aeroelastic ground wind loads model where the maximum dynamic bending moment 
at the base of the vehicle on the launch pad is plotted against wind speed 
in knots. At wind speed of about 57 knots, which corresponds to a Strouhal 
number of about 0.2 based on the diameter of the lower stages, the response 
peaks sharply. At 'a low value of structural damping, £ = 0.01, the time 
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history of response is nearly sinusoidal. This feature implies the existence 
of a self-excited limit -amplitude Oscillation for which the negative aerodynamic- 
damping term and positive structural-damping term are in equilibrium. At higher 
values of structural damping, such as £ = 0.02, the response time history 
changes to a random-amplitude constant -frequency motion which is typical of the 
response of a lightly damped mechanical system to a random forcing function. 

CONCLUSIONS 

Analysis of the results presented from a wind-tunnel study of the mean- 
drag and unsteady lift forces on a stationary and oscillating circular cylinder 
in two-dimensional flow at Reynolds numbers from O.36 million to 18.7 million 
indicates the following conclusions: 

1. The mean-drag coefficient on the stationary cylinder, at Mach numbers 
less than about 0.2, follows the trends established by previous investigators 
and has an approximately constant value of 0-54 for Reynolds numbers between 

4 million and 10 million. 

2. The frequency content of the unsteady lift force on the stationary 

cylinder can be categorized into three regimes dependent upon Reynolds number 
as follows: wide-band random (1.4 x 10^ <j % < 3*5 X 10^), narrow-band random 

( 3*5 X 10^ < Ru < 6.0 x 10^), and quasi-periodie (r n > 6.0 x 10^). 

3. The Strouhal number of the unsteady lift on the stationary cylinder 
in terms of the center frequency of a Strouhal bandwidth follows the trends 
established by previous investigators at Reynolds numbers from 1.4 million to 
8 million, then at hitherto unexplored Reynolds numbers from 8 million to 

17 million is approximately constant at about 0.3* 



4. At Mach numbers less than 0*3> the root -mean- square unsteady lift 
coefficient on the stationary cylinder fluctuates at Reynolds numbers, from 


1.5 million to 8 million, then the range narrows into a single function which 
decreases slowly with higher Reynolds numbers. 

5. A lift force due to cylinder oscillation exists when the cylinder is 
oscillated at or near the aerodynamic Strouhal frequency of the stationary 
cylinder. This lift force increases with increase in amplitude of motion, 
building up to several times the lift on the stationary cylinder. When the 
cylinder is oscillated at frequencies far removed from the aerodynamic Strouhal 
frequency of the stationary cylinder, there is no significant lift due to 
motion. 

6. The unsteady lift due to motion was found to have a destabilizing 
aerodynamic damping component for cylinder motion at frequencies below the 
stationary cylinder vortex- shedding frequency. This component shifts abruptly 
to a stabilizing damping force at frequencies above the vortex- shedding 
frequency. 
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Figure 1. - Comparison of present and previous investigations. 




Figure 2. - Cylinder model installed in the wind tunnel. 
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Figure 6. - Strouhal number of unsteady lift on stationary cylinder. 
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Figure 7- ~ Unsteady lift coefficient on stationary cylinder. 




Figure 8.- PSD f s of unsteady lift without and with cylinder motion. 
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Figure 9 -- Effects of cylinder motion on wide-band RMS lift coefficient. 




